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We study ferromagnetic transition in tliree-dimensional double-exdiange model containing im- 
purities. Tfie influence of both spin fluctuations and impurity potential on conduction electrons is 
described in coherent potential approximation. In the framework of thermodynamic approach we 
construct Landau functional for the system "electrons (in disordered environment) -I- core spins". 
Analyzing the Landau functional we calculate the temperature of ferromagnetic transition Tc and 
paramagnetic susceptibility x- For Tc, we thus extend the result obtained by Furukawa in the 
framework of the Dynamical Mean Field Approximation, with which our result coincides in the 
limit of zero impurity potential. We find, that the alloy disorder, able to produce a gap in density 
of electron states, can substantially decrease Tc with respect to the case of no impurities. We also 
study the general relation between the Coherent Potential Approximation and the Dynamical Mean 
Field Approximation. 
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I. INTRODUCTION 



II. HAMILTONIAN AND CPA EQUATIONS 



The recent rediscovery of colossal magnetoresis- 
tance in doped Mn oxides with perovskite structure 
R-i-ajDajMnOs (R is a rare-earth metal and D is a di- 
valent metal, typically Ba, Sr or Ca) substantially in- 
creased interest in the double-exchange (DE) model . 
Several approaches were used lately to study the thermo- 
dynamic properties of the DE model, including the Dy- 
namical Mean Field Approximation (DMFA) [Q (and ref- 
erences therein; DMFA itself see Green functions de- 



coupling techniques l6 
mean-field approach 



, Schwinger bosons 0| , variational 
I and numerical methods [p|-^d|. 
In all these approaches chemical disorder introduced by 
doping, which is generic for the manganites, has not been 
taken into account. 

Recently we have shown that the concurrent action of 
the chemical and magnetic disorder, is crucial for the the 
description of the density of states and conductivity in 



manganites |12|. In the present paper we consider the fer- 
romagnetic transition in the case of non-zero potential of 
randomly distributed impurities, using the same coherent 
potential approximation (CPA) Jl3[-p^| as in our previ- 
ous paper Q (see also relevant Ref. |1^]). Briefly, the 
effect of impurities on the transition is the reduction of 
Tc as the impurity potential strength increases. In more 
detail, at values of the impurity potential at which the 
electron band is split off into two sub-bands, the mecha- 
nism of ferromagnetic exchange is other than in the case 
of the zero or weak potential. 



We consider the DE model with the inclusion of the 
single-site impurity potential. We apply the quasiclassi- 
cal adiabatic approximation and consider each core spin 
as a static vector of fixed length S {Si = Sn.i, where rij 
is a unit vector). The Hamiltonian of the model H[[\\i\) 
in site representation is 



Hi 



ti-j + (ei - Jrii ■ a) S: 



(1) 



where ti_.j is the electron hopping, is the on-site en- 
ergy, J is the effective exchange coupling between a core 
spin and a conduction electron and a is the vector of the 
Pauli matrices. The hat above the operator reminds that 
it is a 2 X 2 matrix in the spin space (we discard the hat 
when the operator is a scalar matrix in the spin space). 
The Hamiltonian (0) is random due to randomness of a 
core spin configuration [rij] and the randomness of the 
on-site energies e^. 

To handle CPA we present Hamiltonian (|^) as 



H 



Jn,: 



(2) 



(the site independent self-energy S(i?) is to be deter- 
mined later), and construct a perturbation theory with 
respect to random potential V . To do this let us intro- 
duce the exact T-matrix as the solution of the equation 



in which Gq 

G we get 



T^V + VGoT, (3) 
{E — H^)^^ . For the exact Green function 



G = Go -f GqTGq. 



(4) 



1 



The self energy is determined from the requirement 

fc) = Go (5) 

where the angular brackets with the indexes mean av- 
eraging over the configurations of both core spins and 
impurities. In CPA T is considered in a single-site ap- 
proximation, at which Eq. (|^) is reduced to the equation 



y^) =0, (6) 

/ n.e 

where Ti is the solution of the equation 



Here 



f^ = V, + V,g{E)f,. 



giE)^go{E-t), 



(7) 

(8) 
(9) 



No{e) being the bare (i.e. for ei — and J — 0) density 
of states (DOS). Finally Eq. (|) can be transformed to 
an algebraic equation for the 2x2 matrix S: 



(10) 



III. BAND STRUCTURE 



To consider the evolution of the DOS 



N{E) = -luig,, 



(11) 



where gc = g (here Tr means the trace over spin states 
only), with the variation of the impurity concentration 
and potential strength, we exploit the semi-circular (SC) 
bare DOS given at |e| < W {W is half of the bandwidth) 
by 



(12) 



and equal to zero otherwise. For this DOS (which is exact 
on a Caley tree) 



9oiE) = - 



(13) 



Hence we obtain 

± = E-2wg~g~\ (14) 
where w = W'^/S. Thus, Eq. ( p^ ) transforms to 



(E - e + Jna ~ 2wg)~'^) . (15) 



It is convenient to write the locator g in the form 

9 =\{9J + gs^), (16) 

where / is a unity matrix. For the charge locator gc and 
spin locator we obtain the system of equations 



5c = 2. 



gs -2 



E - e - wgc 



{E~e~wgc) -(Jn-wgs)2 



wgs - «/n 



(E-e-wgc) -(Jn-wgs)^ 



(17) 



In the strong Hund coupling limit (J oo) we obtain 
from Eqs. (|l^) two decoupled spin sub-bands. For the 
lower sub-band, after shifting the energy by J we obtain 



9c = 



9s = 



1 



E — e — wgc — wrings 



E - e- wgc 



wn^gs 



(18) 



where gs = (0,0, g^) an axis OZ is directed along the 
average magnetization of core spins m) . 

In fact, details of alloying define how to average over 
the configurations of impurities in Eqs. (0|8|). We use 
for this random substitution model of disorder. That 
is ei = V with the probability x and = with the 
probability 1 — x, where V and x are the impurity po- 
tential and concentration, respectively. As to core spins, 
once CPA is introduced, their configuration probability 
should be determined self-consistently in order to close 
Eqs. ( p7| , p^ . In the Appendix we prove that in the pres- 
ence of an annealed disorder (including dynamical one) 
DMFA ansatz for the disorder configuration probability 
Q keeps the free energy stationary against variations of 
CPA ±. 

For the following two cases, however, only the averag- 
ing over e is left. For a saturated ferromagnetic (FM) 
phase (m = 1), we obtain gs = gc, and closed equation 
for the charge locator is 



9c = 



2wgc 



(19) 



For a paramagnetic (PM) phase (m — 0), we obtain 
gs = 0, and closed equation for the charge locator is 



9i"^ 



1 



E 



(0) 

wgc 
.(0) 



(20) 



It appears that calculation of gc is sufficient for obtain- 
ing Tc and the PM susceptibility x as functions of x and 
V. 



2 



The results of numerical calculation of the DOS at the 
PM state are presented on Fig. 1 (specific x and different 
V). 



V > Vc the gap exists for all < x < 0.5, and a,tV = Vc 
the gap closes exactly for x = 0.5. Anyhow, even if the 
gap is closed there still may exist pseudogap at a strong 
enough potential (see Fig. 2). 

The dependence of the gap A on a; at several V/W is 
shown on on Fig. 3, which summarizes the regularities 
discussed above. 





FIG. 2. The calculated DOS as function of impurity con- 
centration X for V/W = 0.69. 



When comparing Eqs. and (|20|) it is seen that 

the DOS in the FM phase is equal to the DOS in the 
PM state of another model, with increased by a factor 
of V2 bare bandwidth. Using this property and Fig. 1 
we conclude that at an appropriate V/W there may be a 
gap between the conduction band states and the impu- 
rity band states in the PM phase, while the FM DOS is 
gapless. This may explain metal-insulator transition ob- 
served in manganites and magnetic semiconductors [ p^ . 
It is also seen from Fig. 2 that the gap in the DOS exist- 
ing at low concentrations may close at higher concentra- 
tions. The value of the impurity potential Vc — W / \f2 
detaches two types of the gap behavior. KtV <Vc the 
gap opens for some x < 0.5 or does not open at all, at 



It is worth noting, that at F > Vc, A is near to saturate 
after x > 0.25. This may explain the near independence 
of the resistivity activation energy on x in the PM phase 
observed in single crystalline manganites jl^] (the fact 
used by some authors to support the polaron scenarios 
of the transport at T > Tc)- 



IV. LANDAU FUNCTIONAL 

In our approach the calculation of thermodynamic 
properties is based on the analysis of the Landau func- 
tional. We start from the exact partition function of the 
system electrons -I- core spins 

N 



exp 



-/?r!([n,])+/3H. Vn, 



N 

=1 



[]dn„ (21) 



i=l 



where N is the number of core spins, H is magnetic field, 
(3 = 1/T, and ri([nj]) is the grand canonical potential of 
the electron subsystem for a given core spin configura- 
tion [rij]. Our aim is to obtain Landau functional of the 
system. 

Since we use CPA for electronic properties it is consis- 
tent to construct the Landau functional using a mean- 
field approach. Mean-field means that we do not take 
into account large-scale fluctuations of the macroscopic 
magnetization M = -i- J^j : but within CPA we cer- 
tainly take into account microscopic fluctuations of n;. 
In the mean-field approximation all energy levels depend 
only upon M. Hence we may approximately put 



3 



(22) 



V. CALCULATION OF Tc AND x{T) 



where fl{m) is the grand canonical potential calculated 
within CPA at a non-zero m. Then the partition function 
can be written as 



Z 



exp [-/3L(M, H)] dM, 



where 

L(M, H) = n{M) - iVH • M - TS{M); 
the quantity 

SiM) = In / <5 ( M - 1 ^ n, ) n dn, 



(23) 



(24) 



(25) 



is the entropy of the core spin subsystem. So we may 
identify the functional L(M, H) in the exponent of Eq. 
( p3| ) with the Landau functional of the whole system . 

At high temperatures the minimum of i(M, 0) is at 
M = 0. At low temperatures the point M = corre- 
sponds to the maximum of this functional. Let us expand 
L{M, 0) with respect to M 



L(M, 0) = Lo{T) - L2{T)M^ + O (M^), 
L2{T) = fl2{T)-TS2. 



(26) 



Here the coefficients are defined from the the second- 
order expansions of S{M) 



SiM) = So 



where S2 = and of n{m) 

fl{m) = na{T) - Q.2{T)m^ + 



(27) 



(28) 



which is to be constructed. Thus the critical tempera- 
ture Tc, below which the functional has the minimum at 
some M ^ 0, is defined from the equation 



L2[Tc)^^2{Tc)-TcS2 



0. 



(29) 



The magnetic susceptibility in the PM phase is given 
by the equation 



N 



X{T) 



1 



TS2^VL2{T) 



(30) 



The grand canonical potential of the electron subsys- 
tem is given by 



N 



+ C)C 



In 



1 



N(E)dE, (31) 



where N{E) is the DOS given by Eq. ( |ll| ) and fi is the 
chemical potential determined from the equation 



f(E,ti)N{E)dE, 



(32) 



in which f{E,ii) is the Fermi function and n = 1 — a; is 
the number of electrons per site. 



Due to the constancy of n (or x) the calculation of 
5l2(T) requires only the second-order expansion of N{E) 
with respect to m, which is obtained via Eq. ( pT| ) thus 
expanding Eq. ( p^ ) 



(33) 



Substituting the related result for iV^^) (E) into Eqs. 
( ^) and (Hy) , we obtain the following equations 

(34) 



Tc^eiTc,fi{Tc)), 



iV^^^^ T-e{T,fi{T)) 



In these equations /i(T) is determined from Eq. (|3_ 
where N{E) is replaced by its PM value N^^^E), and 



(35) 



e(T,M) 



2w 



In 



Im ■ 



37r/3 J_ 

(9.). [{9'.). + f {(9.) Ad!). -{9!) 




dE, (36) 



where 



9e 



E — e — wg, 



(0) 



(37) 



Using Eq. ( [20| ) and integrating by parts, we obtain 

/oo 
f{E,ti)e{E)dE, (38) 
-00 

where 

?/i I ( n 1 

(39) 



The function 6{E) being integrated with respect to en- 
ergy gives exactly zero. This leads to Tc{x) — Tc{l — x), 
which reflects the particle-hole symmetry of our model 
irrespective of the disorder strength. 

In the case where the DOS is smooth, the inequality 
W ^ Tc allows us to consider electrons as nearly degen- 
erate. In this case both and 0(T, /i) do not depend 
upon T, and Eq. (Q) is just a ready formula for Tc- 
The same is true if there is developed gap. In this case 
fi is near the middle of the gap, so we can substitute 
f{E, /i) by one, provided that the integration in Eq. (pF 
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extends only over the filled sub-band. In both these cases 
x(r) obeys the Curie-Weiss law (see Eq. (|3^)). In the 
first case the ferromagnetic order is mediated mostly by 
mobile holes that is specific for DE. In the second case 
the concentration of the mobile carriers (both holes and 
electrons) is exponentially small. So effective exchange 
between core spins is mostly due to virtual transitions 
of electrons from the lower filled to the upper empty 
sub-band via the gap. This mechanism is an analog of 
super-exchange (SE) acting in the system with electron 
disorder. 

If A ~ r or there is a pseudogap with a strong dip in 
the DOS, the integration in Eqs. (H) and (H) should 
take into account the tails of f{E,iJ,). The exchange in 
such cases is intermediate between DE and SE types. 

For the case of no on-site disorder our Eq. (Q) for 
Tc coincides with Eq. (49) of Ref. obtained in the 
framework of DMFA and calculated numerically versus 
X. In this case we even managed to calculate the integral 
in Eq. (|^ analytically, to get for Tc'- 



WV2 



An 



1 



V2fi/W, 



while Eq. 



takes the form 
1 



TT 



sm 



y + y 



(40) 
(41) 

(42) 



When the disorder is taken into account the integra- 
tion in Eq. (|3^)-(|3^) is done numerically. The results for 
Tc{x) at different V/W are presented on Fig. 4. The 
upper (dashed) curve calculated at = is the same as 
plotted using Eqs. (0)-(|42l) 
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FIG. 4. The calculated Curie temperature as function of x 
at different V/W. 

One notices that the increase of the impurity poten- 
tial leads to progressive decrease of Tc- This decrease 



becomes more substantial at V/W and x able to pro- 
duce the gap in the density of states. This trend marks 
the modification of the ferromagnetic exchange mecha- 
nism which accompanies the electron band splitting. As 
in the case of DE alone Tc still increases from zero to a 
maximum value upon increasing x from a; = to a; = 0.5. 



VI. DISCUSSION 

Our results show that Tq is decreased by the pres- 
ence of impurities with strong enough potential. It may 
be questioned in this connection why low value of Tc 
in manganites would be on account of fluctuations be- 
yond DMFA in the pure DE model rather than due to 
the alloy disorder. The present study reveals an inter- 
esting, though hardly experimentally detectable, feature 
- the alternation of ferromagnetic exchange mechanism 
from DE to SE like as the impurity band splits off the 
conduction band. 
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APPENDIX A: DYNAMICAL MEAN FIELD 
APPROXIMATION FOR GENERAL STRONGLY 
CORRELATED SYSTEM 

In most interesting cases the partition function of a 
strongly correlated electron system may be represented 
as functional integral 



Z = 



I Z [v, h] H [n,] Ph [h.] Vv, (r) Ph, (r) (Al) 



in which the integrand Z [w, h] is the partition function 
of non-interacting electrons system posed in imaginary- 
time T and site i dependent potential Vi (r) and mag- 
netic hi (r) fields distributed independently on different 
sites with on-site probability densities Py [v] and Ph [h], 
respectively. The free-electron partition function is rep- 
resented by the integral 



Z[v,h] 



exp 




d_ 



+ H-n{T) 



(A2) 



over anticommuting fields ip = {"(Aia} and its conjugate 
"0, where the imaginary-time t dependent Hamiltonian 
Ti, (r) is defined by the matrix elements 



5 



Vq/3 (i; t) = Vi (r) Sa,f3 - hi (r) • aa/3 (A3) 

(the electron charge and magneton are absorbed in the 
corr esponding fields). The Grassmanian integral in Eq. 
( [A^ ) is calculated in the familiar form 

Z[v,h] =Detg-^ =exp{-S[v,h]} (A4) 
S[v,h] = log Det^; = Tr log G (A5) 

where Q^^ = + /j, — Ti,\ the matrix element of Q with 
respect to site, spin and imaginary time, Gia.jfs {t,t'), 
is the fermionic Green-Matsubara kernel satisfying the 
equation 

d 



Q:;: + ^^j Sia.jfi (t, t') 

■ ^ Ti-iaM-i ) Sk-f,]f3 (t, t') = Iia,jf3 [t, t') (A6) 
A;, 7 



where 



Iia,j0 (r, r') = S,jSa,f35 (r - t') 



(A7) 



is the unity matrix. Thus we obtain for the partition 
function 

Z= flip, [v,] Ph [h.] exp {-S [v, h]}Vv, (r) Vh, (r) 

(A8) 



and for the number of electrons 
N, 



(A9) 



(which is the equation for chemical potential ^) where 
the average of every functional $ = (f> [u, h] is defined by 



exp {-S [v, h]} $ [i;, h] Vv^ (r) Vh, (r) 



(AlO) 



Assume that, in addition to 'annealed' disorder described 
by Vi{T) and hi^r), there exists a 'quenched' disorder 
described by independently fluctuating static on-site en- 
ergies e^. Then the fields Vi (r) in the Hamiltonian are 
to be replaced by Vi (r) -I- Ci and the observables should 
be additionally averaged over given distribution of the 
quenched disorder. For example, free energy is given by 



F 



-^log^^ +AiiV, = i(-logZ-HMTr(g)„j^ 



(All) 



Gonsider an effective medium described by the 
bare hopping Hamiltonian t plus some self-energy 



given by translationally invariant and stationary matrix 
{i — j,T — t'). The Green-Matsubara kernel for such 
a medium is also translationally invariant and stationary: 
GiajfS (t, t') — Gap {i — i,T — t'), and can be expressed 
via Fourier transform 



G (k, 



/i — i (k) — S (k, iuj„ 



(A12) 



where ujm is the Matsubara frequency, Eq/j (k, iiOm) is 
Fourier transform of Sa/3 {i — j,T — t') and the quanti- 
ties in both sides of Eq. ( A12) are matrices in spin space. 
From the very meaning of effective medium it is to be re- 
quired 



(^>,,.h 



(A13) 



that implicitly defines the self-energy. Eq. (A13) can be 
transformed to another form using T-Matrix ansatz for 
the exact Green-Matsubara kernel. We have 



G = G + GTG 



where T is the exact T-Matrix 



T={I-VG) ' V 



(A14) 



(A15) 



where V is the effective interaction with the matrix ele- 
ments 

Via j73 (t, t') = Via,j0 (t) (5 (t - t') - Sia.jp (t, t') (A16) 

(Note that the effe ctive interaction unlike true one is non- 
local). Then Eq. ( A13 ) is exactly equivalent to following 
equation 



(^>.,h 







(A17) 



which expresses the requirement of zero scattering on av- 
erage. It can easily be understood that the complexity 
of solving Eq. ( A17) to determine (z — j, t — r') or 
(k,iujm) is the same as that of calculating G at ev- 
ery disorder configuration and then performing the aver- 
age. To develop a working approximation let us apply the 
idea of Goherent Potential Approximation (GPA) (used 
exceptionally for quenched disorder) to the present case. 
Within GPA one firstly puts 



{i ~ j,T - t') = dij'SaP {t - t') 



(A18) 



(so V becomes also local in space) and secondly applies 
the zero-order locality approximation 



Gal3 {i - j,T- t') ~ gafS (t - t') (5,; j 
9a0 (t) = Gap (0,t) 



(A19) 
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to the term VG in Eq. (A15). This leads to approximate 
space locahty of the transfer matrix as foUows 



%aji3 (t, r') w Tap {i; r, r') 6^, 



(A20) 



where the elements of the local transfer matrix Tj satisfy 
the integral equation 



/3 



[Vaj {i; t)S{t- v) 



-Sq^ (t - v)] g^s {v - v') Tsp (i; v, v') dvdv' 

= (i; r) <5 (r - r') -S„;3 (r-r) (A21) 



The relaxed form of Eq. (A17), which is essentially CPA 
equation, reads 



T,) ) =0 



We use symbolic solution of Eq. ( A21| ) 



T,; = 



(A22) 



(A23) 



where /, Vi, E and g are the matrices with the ele- 
ments SapS {t-T^, Vap {i;T)d (j-T^, Sq^ (^^'^') 

and gap ^ ^ respectively. Hence one may rearrange 
Eq. ( |A22[ ) to the form 



(A24) 



There arises a question of how to consistently perform av- 
erage over V and h in Eqs.( [A2^ ),( A24 ), which may sound 
strange because the answer seems ready. Indeed, upon 
the approximation given by Eq. (A20) and the exact 
decomposition given by Eq. (A14) one gets the approx- 
imate expression for Q which would be used to calculate 
the approximate entropy functional due to Eq. (|A5|): 



S* [w, h] ~ log Detg + log Det (X + TQ) 



(A25) 



However this expression is still too cumbersome to be 
used for averaging since Q has intersite matrix elements, 
neglected when solving for T (it concerns only the second 
term because the first doesn't depend on random fields). 
But if we again apply the approximation of Eq. ( Al£ ), 
this time to the term TQ, and make use of Eqs. (A23), 
( A24| ) we obtain just site-additive entropy functional 



5[i;,h] -Trlog^ + ^s, K,h,], 

i 

Si [vi, hj] = log det {j: + fig^ 



(A26) 



where 'tr' and 'det' mean corresponding operations only 
with respect to spin and r variables. The admission of 



S [w, h] in the form of Eq. ( |A26[ ) is just Dynamical Mean 
Field Approximation (DMFA) ansatz. Let us show that 
DMFA is the best in variational sense, consistent with 
CPA, approximation, namely that it ensures the station- 
arity of the approximate free energy 

/3"^^/log / exp{-s, [u,;,h,]} 

i ^ 

Pv K] Ph [h.] Vv, (r) Vh, (r))^ (A27) 

with respect to any variations of the local self-energy. To 
prove this let us first note that variation of the first two 
terms in brackets on account of varying /i cancels out. So 
we obtain upon variation of the matrix E 



g-(((g-' + ±-V, 



(A28) 



where Sg~^is the variation of the matrix g^^ due to the 
variation of E. Now it is seen that for the stationarity, 
that is SFa = 0, it is necessary and sufficient satisfying 



Eq. (A24) or equivalent Eq. (A22) 



APPENDIX B: DOUBLE EXCHANGE MODEL 
WITH CLASSICAL CORE SPINS 

In the case of DE Model considered we have Vi (t) = 
and static magnetic fields 



So 



hi (r) = Jui 



PM=Hnt-l) 



(Bl) 



(B2) 



In present case we can use Fourier transform with respect 
to T. We obtain the basic equations 

/ r cxp{~..[e,n]} 



/ exp {— s [e, n]}c?n 



where 



s [e, n] = ^ logdet |/ + E (iw„) — e + Jn • ct g (iujn)^ 

(B4) 



and 



g (iuJm) ^ No (e) iw„ -I- ^ - e - E (iuJm) de. (B5) 
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